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4.1 Introduction 


Low-level warm marine boundary layer (MBL) clouds cover large regions of Earth’s 
surface. They have a significant role in Earth’s radiative energy balance (Klein and 
Hartmann, 1993) and hydrological cycle. Despite the fundamental role of low-level 
warm water clouds in climate, our understanding of these clouds is still limited. In 
particular, connections between their properties (e.g. cloud fraction, cloud water 
path, and cloud droplet size) and environmental factors such as aerosol loading and 
meteorological conditions continue to be uncertain or unknown. Modeling these 
clouds in climate models remains a challenging problem. As a result, the influence 
of aerosols on these clouds in the past and future, and the potential impacts of these 
clouds on global warming remain open questions leading to substantial uncertainty 
in climate projections. To improve our understanding of these clouds, we need 
continuous observations of cloud properties on both a global scale and over a long 
enough timescale for climate studies. At present, satellite-based remote sensing is 
the only means of providing such observations. 

The cloud droplet effective radius (r e ) is one of the most important cloud para- 
meters that are routinely monitored from space. The r e is defined as (Hansen and 
Travis, 1974) 


J 0 °° r 3 n(r, z ) dr 
J 0 °° r 2 n(r, z) dr ’ 


(4.1) 


where n(r, z) is the cloud droplet size distribution (DSD) at altitude z in cloud. The 
cloud droplet effective radius determines the optical thickness of cloud for a given 
amount of water (Twomey, 1974). It also has a significant influence on precipitation 
formation processes (Lebsock et ah, 2008; Kubar et ah, 2009). Therefore, it is 
a key microphysical parameter required to estimate radiative effects of clouds, 
study aerosol-cloud-precipitation interactions, and validate cloud parameterization 
in global climate models. 

Many remote sensing methods exist to infer cloud r e from various types of 
satellite instruments (e.g. Prabhakara et ah, 1988; Nakajima and King, 1990; Austin 
et ah, 2009). Of particular interest in this review are cloud r e retrievals from MODIS 
(or Moderate Resolution Imaging Spectroradiometer). MODIS is a key instrument 
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aboard the Terra and Aqua satellites. Terra’s orbit around Earth is timed so that 
it passes from north to south across the equator in the morning, while Aqua passes 
south to north over the equator in the afternoon. Terra and Aqua MODIS are 
viewing the entire Earth’s surface every one to two days, acquiring data in 36 
spectral bands. A variety of cloud parameters, from cloud fraction to cloud top 
height to cloud optical thickness (r) and r e , can be inferred from MODIS multi- 
spectral observations (Platnick et al., 2003). The operational MODIS cloud retrieval 
algorithm, as described in detail in the next section, is based on the so-called 
bi-spectral solar reflective method (referred to as bi-spectral method hereafter), 
which utilizes cloud reflectance measurements from two spectral bands to retrieve 
cloud r and r e , simultaneously (Nakajima and King, 1990). One measurement is 
usually made in the visible or near-infrared spectral region (e.g. 0.86 pm), where 
water absorption is negligible and therefore cloud reflection is mainly determined 
by r, and the other in the shortwave infrared (SWIR) (e.g. 2.1 pm or 3.7 pm), 
where water is significantly absorptive and cloud reflectance primarily decreases 
with increasing cloud droplet size. The r and r e retrievals based on the bi-spectral 
method are widely used for validating climate models (Kay et ah, 2012; Pincus et 
ah, 2012), studying aerosol-cloud interactions (Quaas and Boucher, 2005; Quaas 
et ah, 2009) and facilitating other cloud remote sensing techniques (Lebsock and 
L’Ecuyer, 2011). 

Unlike cloud optical properties (such as r) that can be spectrally dependent, 
cloud r e , a physical cloud parameter, should be independent of the observation 
method or what spectral band is used. For example, the MODIS instrument has 
three SWIR bands, centered at 1.6 pm, 2.1 pm, and 3.7 pm, respectively. In the 
MODIS operational cloud product (MOD06) Collection 5 (C5) processing algo- 
rithm, the combination of the 0.86 and 2.1 pm bands are used for r and r e retrievals 
over open ocean (Platnick et ah, 2003). Hereafter, we will refer to the r e retrieval 
based on the 2. 1-pin band observation as r e , 2 .i- Besides r e , 2 . 1 , the MODIS also 
provides two other r e retrievals, one based on the 1 . 6 -pm band and the other based 
on the 3.7-pm band observations (hereafter referred to as r e p . 6 and r ej 3 . 7 ). One 
might expect these three r e retrievals to be in close agreement. However, several 
studies have found substantial differences between them (Nakajima et ah, 2010a; 
Seethala and Horvath, 2010; Zhang and Platnick, 2011). In particular, it is found 
that: 

MODIS r ei 3.7 retrievals for marine warm clouds are generally smaller than r e ^.i 
(and r e ,i. 6 ); 

Geographically, the differences between r e ^ 7 and r e 2 .i based on MODIS Level 3 
monthly mean product show strong dependence on cloud regimes, small varia- 
tions (close to zero) over the costal stratocumulus regimes, and large variations 
over the cumulus cloud regimes; 

At pixel level, the difference between r e , 3.7 and r 6 i 2.1 correlates with cloud r, 
r e , and the degree of sub-pixel inhomogeneity. 

A detailed analysis of the differences between r e , 3.7 and 7 ^ 2. 1 MODIS retrievals will 
be given in section 4.3. r e , 1.6 is not considered in this study, mainly because the 
1 . 6 -pm band on Aqua MODIS has nonfunctional or noisy detectors that lead to the 
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striping issue (Wang et al., 2006). As a result, it is difficult to make pixel-to-pixel 
comparisons between r e> i.g and other effective radius retrievals. 

The spectral dependence of MODIS r e raises many questions. The most impor- 
tant one, for all the users of MODIS products, is whether the difference between 
r e , 3.7 and r e ,2.i is an artifact due to, say, algorithm issues and inherent limitations 
of the retrieval method, or something meaningful that contains information about 
the cloud. Several lines of evidence suggest that it is unlikely to be due to technical 
issues (e.g. code bugs, ancillary data issues, etc.) of the operational MODIS cloud 
retrieval algorithm. First, several studies based on independent retrieval algorithms 
also found substantial differences between r ej 3.7 and r ei 2.1 that are similar to the 
operational MODIS cloud product (Nakajima et ah, 2010a; Minnis et ah, 2011). 
Second, both theoretical and numerical studies indicate that r e> 3.7 and r e ^.i are 
expected to differ significantly under certain circumstances owing to their difference 
sensitivities to, for example, cloud vertical structure (Platnick, 2000) , 3D radiative 
effects (Zhang and Platnick, 2011; Zhang et ah, 2012), and the presence of drizzle 
drops in the cloud (Nakajima and King, 1990; Zhang, 2013). 

The question is then: what has caused the differences between r 6 ) 3.7 and r 6 i 2 . 1 ? 
A number of recent studies have attempted to address this question from dif- 
ferent perspectives (Platnick, 2000; Nakajima et ah, 2010a, 2010b; Seethala and 
Horvath, 2010; Zhang et ah, 2010, 2012; Zinner et ah, 2010; Zhang and Platnick, 
2011; Zhang, 2013). Several mechanisms have been proposed to explain the above- 
mentioned spectral dependence of MODIS r e retrieval, which can be divided into 
two categories based on their underlying physics and consequent implications. In 
the first category are those mechanisms related to the fact that r e , 3.7 and r 6 j 2.1 have 
different sensitivities to cloud vertical structure and the presence of large drizzle 
drops in warm water clouds, with implications that the spectral difference actually 
carries useful information about the cloud that can be used for remote sensing 
and model validation. In the other category are those mechanisms related to the 
inherent limitations of the bi-spectral method, such as, the lack of consideration of 
3D radiative effects and sub-pixel inhomogeneity in the retrieval method. In such 
cases, the spectral difference implies significant retrieval uncertainties should be 
cautioned when using MODIS r e retrieval products. Overall, these studies reveal 
that the spectral difference of MODIS r e retrievals is a complicated issue that has 
complex causes. 

This chapter provides an overview of the current understanding of the spectral 
dependence of MODIS r e retrievals. In section 4.2, we briefly revisit the opera- 
tional MODIS cloud r e retrieval algorithm to set the stage for later discussion. In 
section 4.3, the MODIS r e> 3.7 and r 6 ) 2.1 retrievals for MBL clouds are compared to 
illustrate the spectral dependence of MODIS cloud r e retrieval. Section 4.4 intro- 
duces several potential mechanisms that may cause or contribute to the spectral 
dependence of MODIS r e retrieval. The relative role of these mechanisms in dif- 
ferent cloud regimes is discussed in section 4.5. Current outstanding issues and an 
overview for future work are given in section 4.6. 
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4.2 Operational MODIS r e retrieval algorithm 


In order to set the stage for later discussion in this section, we briefly describe the 
operational MODIS r e and r retrieval algorithm based on the bi-spectral method. 
The retrieval makes use of a pair of measurements of cloud reflectivity, one from a 
visible or near-infrared (VNIR) MODIS band (0. 86-prn band over ocean) and the 
other from a shortwave infrared (SWIR) band (e.g. 2.1 pm or 3.7 pm band) (Naka- 
jima and King, 1990; Platnick et ah, 2003). The visible band measurement provides 
the information for r retrieval because water absorption in the VNIR region is al- 
most negligible and, as a result, cloud reflectivity is mainly determined by r. The 
SWIR measurement provides the information for r e retrieval because significant 
water absorption in SWIR makes cloud reflectance decrease with increasing cloud 
r e . In the operational algorithm, this method is implemented by using the so-called 
look-up-table (LUT), as shown in Fig. 4.1. The LUT contains pre-computed cloud 
reflectivities at visible and SWIR bands for various combinations of r e and r under 
different Sun-satellite viewing geometries and surface reflectances. As illustrated in 
Fig. 4.1, in practice. r e and r are retrieved by projecting the observed reflectivities 
onto the LUT. Once r e and r are retrieved, the liquid water path (LWP) of the 
cloud can be easily derived from the equation LWP — | prr e assuming that the 
cloud layer is vertically homogenous and the extinction coefficient Q e of a cloud 
droplet is about 2. The MODIS cloud r e and r retrievals have been shown to agree 
reasonably well with other satellite products (e.g. Zhang et ah, 2009; Minnis et ah, 
2011; Pincus et ah, 2012; Stubenrauch et ah, 2012; Walther and Heidinger, 2012) 
and in situ measurements (e.g. Painemal and Zuidema, 2011; King et ah, 2012). 


a) R (0.86 pm) & R (2.1 pm) LUT 



b) R (0.86 pm) & R (3.7 pm) LUT 



Fig. 4.1. Examples of the look-up-table of cloud bi-directional reflection function as 
functions of cloud optical thickness and effective radius, based on the combination of (a) 
0.86 and 2.1-pm bands, and (b) 0.86 and 3.7-pm bands. 
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It is important to note that the bi-spectral method is based on several important 
assumptions about the cloud (or the cloudy pixel): 


1 . 

2 . 

3. 


Within the pixel, the cloud is vertically homogeneous, (referred to as the vertical 
homogenous assumption). 

Within the pixel, the cloud is horizontally homogenous. The cloud reflectance 
of the pixel of interest is independent of the properties of surrounding pixels, 
(referred to as the plane-parallel and homogeneous cloud assumption). 

Within the pixel, the sizes of cloud droplets follow certain analytical distri- 
butions, namely the following monomodal Gamma distribution (King et ah, 
1998): 


1 — 3v e 

n(r ) = Nr v « exp 



(4.2) 


where n(r ) is DSD, N is a constant, and v e is the effective variance assumed 
to be 0.1 in the operational MODIS algorithm (monomodal Gamma DSD as- 
sumption hereafter) . 


These assumptions are necessary because, at most, two independent pieces of in- 
formation can be extracted from a pair of cloud reflectivities. Since the information 
content has been used for r e and r retrievals, other aspects of cloud would have 
to be assumed. These assumptions are justified for some clouds, particularly, non- 
precipitating stratocumulus clouds (Martin et ah, 1994; Di Girolamo et ah, 2010). 
However, for other clouds, they can be problematic. For example, it is known that 
trade wind cumulus clouds can be far from plane-parallel (Liang et ah, 2009; Di 
Girolamo et ah, 2010). The vertical homogeneous assumption and the monomodal 
Gamma DSD often break down when precipitation begins to form within the cloud. 
The warm rain processes, such as collision-coalescence, could broaden the DSD and 
sometimes give rise to a second mode, the so-called drizzle or rain mode, leading to 
a bi-modal DSD (Berry, 1967; Berry and Reinhardt, 1974; Pruppacher and Klett, 
1997). In addition, the collision-coalescence processes make drizzle drops grow big- 
ger as they fall from cloud top towards cloud base, inducing vertical structures 
within the cloud (Pruppacher and Klett, 1997). When these conditions occur and 
the above-mentioned assumptions break down, MODIS r e and r retrievals face 
substantial uncertainties. It will be shown later in section 4.4 that the mechanisms 
that cause significant differences between r e> 3.7 and r e> 2.1 are more or less connected 
to breakdown of these fundamental assumptions about cloud made in MODIS re- 
trieval. 


4.3 Spectral dependence of MODIS r e retrievals for 
MBL clouds 

As stated from the beginning, we focus only on the warm liquid-phase clouds over 
ocean. Since such clouds reside mostly in the MBL, we will refer to them as MBL 
clouds for simplicity. We attempt to identify a MODIS pixel as MBL cloud based 
on the following criteria: (i) the pixel is labeled as ‘confident cloudy’ by the 1 -km 
MODIS cloud mask product (MOD35); (ii) over ocean; (iii) labeled as ‘liquid water’ 
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by the MODIS 1-km ‘Cloud Phase Optical Properties’ data set within MOD06; 
(iv) cloud top temperature warmer than 273K. These conditions are expected to 
remove most of the situations that may complicate the analyses, such as thin cirrus 
overlapping lower clouds. 

4 . 3.1 Geographical pattern 

The monthly mean r e 2.i and r e> 3.7 for MBL clouds based on May 2007 MODIS/ Aqua 
data are shown in Fig. 4.2a and 4.2b, and their difference Ar e ^. 7 - 2.1 = ?"e, 3.7 — r e, 2.1 
in Fig. 4.2c respectively. Several well-known coastal stratocumulus regions, such as 
off the coasts of California, Peru, and Namibia/ Angola, are clearly seen in the figure. 
One can also see from Fig. 4.2 that the transition from the coastal stratocumulus 
cloud regimes to the offshore cumulus cloud regimes is quite sharp. Accompanying 
this transition, r e ^.i increases substantially from 8 ~ 10 /mi near the coast to as 
large as 20 ~ 25 /mi far offshore. However, r e 3.7 is significantly smaller than its 
counterparts. It is easily seen that Z\r 6j 3 . 7 - 2.1 shows an obvious dependence on 
cloud regime. For example, over coastal stratocumulus cloud regions, Ar e , 3 . 7 - 2 .! 
is close to zero, or even slightly positive. However, over the broken cumulus cloud 
regions, where water cloud fraction is small, r e 3.7 is seen to be smaller than r e ^.i 
by as much as 5-10 //in on average. 


4 . 3.2 Correlation with key cloud parameters 

In the rest of this section, we explore the correlations between MODIS r e retrieval 
differences and several key factors. In doing so, we attempt to identify regimes 
where the r e retrieval differences can be attributed more to cloud physics, such 
as drizzle, than to retrieval uncertainties and artifacts caused by, for example, 3D 
radiative effects. 

Figure 4.3a shows the joint histogram between r 6j 3.7 and r e; 2.1 based on about 
1.5 billion Level 2 marine water cloud pixels collected during May 2007 by Aqua 
MODIS from 60S to 60N. Evidently, the density of points is highest along the one- 
to-one line, attesting that pixel-level r 6j 3.7 and r e , 2.1 modes agree reasonably well. 
It is interesting to see that the bias between r 6i 3.7 and r e , 2 .i is quite small when r e is 
smaller than about 12 ~ 13 /mi. However, when r e is larger than about 15 /mi, the 
histogram distribution is clearly weighted toward the r ej 2 .i side and the deviation 
from the one-to-one line increases with increasing r ei 2 .i- Figure 4.3b presents the 
same story, but from a different perspective. It shows the joint histogram between 
^ r e, 3 . 7 — 2.1 and r ei 2 .i- Note that, in order to reduce the data sampling impact on 
the shape of the joint histogram, the histogram of Ar ej 3. 7 - 2.1 at r e , 2 .i bin has been 
normalized with respect to its peak value at that bin. Therefore, the red color in 
Fig. 4.3b corresponds to the most frequent Ar St 3. 7 - 2.1 at a given r 6i 2.1 bin. The 
gray dotted line in Fig. 4.3b shows the PDF of r e , 2 .i derived from one month of 
Level 2 data. Interestingly, the bin-normalized histogram of Ar et 3. 7 - 2.1 vs. r ei 2 .i 
not only shows an increasing spread in Ar £t 3 . 7 - 2.1 with increasing r e , 2 . 1 , but also 
a clear systematic transition, in terms of the most likely observed Ar e , 3 . 7 - 2.1 for a 
given r ej 2 .i (i.e. the red area), from near-zero values when r e: 2.1 < 15 /mi to larger 
negative values when r ei 2 .i > 15 /mi. The fact that this threshold-like behavior 
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(a) r e 2A [jim] 



Fig. 4.2. Monthly mean MODIS (a) r e , 2 .i, (b) r ej 3 . 7 , and (c) Ar e ^.T-i.\ retrievals for 
warm (cloud top temperature larger than 273 K) liquid-phase clouds over ocean. 
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(a) (b) 



5 10 15 20 25 30 5 10 15 20 25 30 

A, 2.1 [Hm] r e2A [jLim] 

Fig. 4.3. The joint histogram of pixel level (a) r e , 2 .i, vs. r e , 3 . 7 , (b) Ar e 3.7 — 2.1 vs. r e , 2 .i, 
which is bin-normalized with respect to the PDF of r e , 2 .i (gray dashed line). 

takes place at r e ~ 15 /im is particularly interesting because r e ~ 15 gm has 
been suggested to be the threshold for the collision-coalescence process to take 
place in marine water clouds (Gerber, 1996). However, further studies are needed 
to determine whether this is simply a coincidence or due to more fundamental 
physical reasons. 

Figure 4.4 shows the bin- normalized joint histograms between Ar e ^_ 7 - 2.1 and 
cloud optical thickness (r). One can easily note that, when clouds are optically thin 
(e.g. r < 5), Ar e , 3 . 7 - 2.1 varies quite remarkably from —15 /im to 10 gm. However, 
when the cloud becomes sufficiently thick (r > 5), the statistics of Ar &) 3 . 7 - 2.1 
become stable and show little dependence on r. The large variability of Ar ey 3 . 7 - 2.1 
for thin clouds in Fig. 4.4 is likely a result of the combined effects of random retrieval 
uncertainties (see section 4.4.1), 3D radiative effects (section 4.4.5), and the so- 
called plane-parallel r e retrieval bias (section 4.4.4). For thin clouds, the signal is 
comparable or smaller than, the noise caused by instrument uncertainties, ancillary 
data uncertainties, and discretization and interpolation of the LUT. As a result, 
the uncertainty associated with the MODIS r e retrievals for thin clouds is large. 
In addition, thin and broken clouds often have significant horizontal heterogeneity, 
providing favorable conditions for 3D radiative effects and plane-parallel r e retrieval 
bias. Caution must therefore be taken when interpreting the meaning of Ar e ^. 7 - 2.1 
for clouds with r < 5 because many factors other than cloud physics, such as 
retrievals errors and artifacts, all play a significant role in this regime. 

When analyzing the potential connection between cloud horizontal inhomogene- 
ity and r e retrieval failure, we will use a so-called sub-pixel inhomogeneity index 
(H a ) from MODIS. It is defined as (Liang et al., 2009) 

stdev[Aj(0.86 /.im, 250 nr)] 

" nrean[Aj(0.86 /.im, 250 nr)] 


(4.3) 
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0.1 1.0 10.0 100.0 

x [unitless] 


Fig. 4.4. Joint histogram of Ar e , 3 . 7 - 2. 1 and cloud optical thickness (r), bin-normalized 
by the PDF of r. 


where stdev[J?j(0.86 pm, 250 m)] and mcan[i?, 7 ;(0.86 pm, 250 m)] indicate the 
standard deviation and mean of the measured reflectances, respectively, for the 
principle 16 250-m-resolution sub-pixels within the 1-km MODIS pixel retrieval 
footprint. Thus, H a has a spatial resolution (i.e. 1 km) consistent with the cloud 
property retrieval and increases with pixel inhomogeneity. Recent studies found 
that the H a index derived from high-resolution (250-m) MODIS cloud reflectance 
measurement provides an objective and quantitative measurement of the horizontal 
inhomogeneity of 1-km MODIS pixel (Di Girolamo et ah, 2010; Zhang and Platnick, 
2011; Zhang et ah, 2012). The H a index will be reported in the MOD06 product 
in the coming Collection 6. The dependence of Ar e ^. 7 - 2.1 on cloud horizontal in- 
homogeneity is shown in Fig. 4.5. Figure 4.5a and 4.5b show the bin- normalized 
joint histograms of the sub-pixel cloud inhomogeneity (H a ) defined vs. r e , 2.1 and 
r e; 3 . 7 , respectively. Figure 4.5c shows the bin-normalized joint histograms of the 
H c , vs. Ar £} 3 . 7 - 2 .i- Optically thin clouds (r < 5) are excluded from this figure for 
the above-mentioned reason, but results are similar if we include thin pixels (not 
shown). The most compelling feature in Fig. 4.5a is the sharp transition of r e , 2.1 
at H a around 0.3 ~ 0.5. When H a is smaller than 0.3, the most likely r e , 2.1 (he. 
red area in Fig. 4.5a) for a given H a stays relatively constant, within 10 ~ 15 pm. 
However, when H a exceeds about 0.3, the most likely r e , 2 .i value increases dra- 
matically with Hr,. Interestingly, this is not the case in Fig. 4.5b, where the most 
likely value of r e , 3.7 shows only weak dependence on the sub-pixel inhomogeneity. 
It is therefore not surprising to see in Fig. 4.5c the most likely values of Ar e ^. 7 - 2.1 
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Fig. 4.5. Joint histograms of (a) r e , 2.1 vs. sub-pixel inhomogeneity index (H a ), (b) r e , 3.7 
vs. H^, and (c) Ar e , 3. 7-2.1 vs. H a . All histograms in this figure are bin-normalized with 
respect to the PDF of Note that thin clouds with r < 5 have been excluded in this 
figure, but results are similar if we include thin pixels (not shown). 


shifting from near-zero to the negative side when H a exceeds about 0.3. The po- 
tential reasons for the dependence of Ar e ^ 3. 7-2.1 on cloud inhomogeneity H a will 
be discussed later in sections 4.4.4 and 4.4.5. 
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In addition to r e , t, and H a , we have also investigated the dependence of 
Ar e> 3 . 7 _ 2 .i on other factors such as cloud top temperature, solar zenith angle, 
satellite viewing angle, scattering angle, surface reflectance, etc., none of which 
shows an impact on Z\r e , 3 . 7 - 2.1 as dramatic and clear as r e , r, and H a . 


4.4 Potential reasons for the spectral difference 

The fact that Ar e> 3. 7-2.1 is correlated with several key cloud parameters, including 
r e , r, and H a , indicates that the spectral dependence of the MODIS r e retrieval is 
a complex issue likely caused by multiple mechanisms. Indeed, several hypotheses 
have been proposed to explain the causes of Ar e , 3 . 7 - 2 .! and its behaviors described 
in the previous section. This section provides an overview of these hypotheses. It is 
helpful to begin with a classification. The existing hypotheses can be grouped into 
two categories based on their underlying mechanisms and consequent implications. 
In one category are those related to the inherent limitations of the bi-spectral 
method, such as retrieval uncertainties (section 4.4.1), the sub-pixel inhomogeneity 
(section 4.4.4), and lack of consideration of 3D radiative effects (section 4.4.5) in 
the retrieval method. In such cases, Ar e ^. 7-2.1 carries little, if any, information 
about the microphysical property of the pixel, but is rather a retrieval artifact 
indicating significant uncertainties in either r e; 2.1 or r e> 3.7, or both. In the other 
category are those related to the fact that r e< 3.7 and r e< 2.1 have different sensitivities 
to cloud vertical structure (section 4.4.2) and the presence of large drizzle drops 
(section 4.4.3) in warm water clouds. The implication of this is that the spectral 
r e difference actually contains useful information about the cloud that can be used 
and actually has been used for remote sensing and model validation. 


4.4.1 Random error 

As shown in Fig. 4.4, Ar e> 3. 7-2.1 for clouds with r < 5 varies more widely than 
that for thicker clouds. The large variation of Ar et 3. 7-2.1 for thin clouds is likely 
caused by random retrieval uncertainties. There are many sources of uncertainty 
in the MODIS retrieval, such as instrument uncertainties, ancillary data uncertain- 
ties, and discretization and interpolation of the LUT. These sources are usually 
uncorrelated, leading to random errors in the observed cloud reflectances. When 
the cloud is thin, the signal from the cloud is comparable to, or even smaller than, 
the retrieval uncertainties, which can result in large errors in r e , 2.1 and r ei 3.7, and 
a highly variable Z\r e> 3. 7-2.1- An example to illustrate the impact of random un- 
certainties on MODIS t and r e retrievals is given in Fig. 4.6. In the example, one 
point in the MODIS operational LUT, with r = 4.1 and r e — 16 /im, (i.e. the cen- 
ter of the cross in Fig. 4.1), is chosen for the purpose of illustration. Uncorrelated 
random errors, which are assumed to follow the normal distribution with standard 
deviation a — 10%, are added to the reflectances of this point at all three bands to 
mimic the above-mentioned retrieval uncertainties. The magnitude (one standard 
deviation) of error is indicated by the size of the crosses in Fig. 4.1. One million 
such samples are generated and then used to obtained r and r e retrieval samples. 
The probability density functions (PDF) of r retrievals are shown in Fig. 4.6a. 
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The two r retrievals, one based on 0.86-//m and 2.1-/im band LUT and the other 
based 0.86-//m and 3.7-//ni band LUT, both vary closely around the true value 
of r = 4.1 and they agree with each other very well. The PDFs of r e> 2 .i, r e; 3 . 7 , 
and Z\r e> 3 . 7 _ 2 .i are shown in Fig. 4.6b. Interestingly, r ei 2 .i has a narrower PDF 
than r ei 3 . 7 , even though the same magnitude of error (i.e. 10 % in reflectance) is 
given to the two bands. This difference is mainly due to the fact that the 0.86- 
/im and 2.1-/im band LUT (LUT 0 . 86 & 2.1 hereafter) is generally more condensed 
in the r e direction than the 0.86-//m and 3.7-qm band LUT (LUT 0 . 86 & 3 . 7 ) (i.e. 
|<91n7?(2.1 /im)/dr e \ < |<91n7?(3.7 /xm)/<9r e |). As a result, the same magnitude of 
reflectance error (e.g. 10%) leads to different r e error. The fact that LUT 0 . 86 & 2.1 is 
generally less orthogonal than LUT 0 . 86 & 3.7 also plays a role. The orthogonality of 
LUT is determined by the extent to which the cloud reflectance in the SWIR band 
is independent from that in the 0.86-/im band. Because is generally more orthogo- 
nal than (i.e. \dR(2.1 /rm)/<97?(0.86 /iin) | > \dR(3.7 /j,m)/dR(0.86 /xm)|), error in 
the 0.86-/im band has less impact on r e , 3.7 than r e , 2 .i, which contributes to the 
PDF difference between r ej 2 .i and r e , 3.7 in Fig. 4.6b. This effect of LUT orthogo- 
nality is demonstrated in Fig. 4.6c and 4.6d, which shows the PDFs of r retrievals 
(Fig. 4.6c), r e , 2 .i, A:, 3 . 7 ; and Ar e> 3 . 7 - 2.1 (Fig. 4.6d) when we only add random errors 
to 77(0.86 /im) but not to 77(2.1 /im) and 77(3.7 / 1 m). Due to the LUT orthogonality 
difference, the error, which is solely in the 77(0.86 /im) , causes only a small error 
in r e , 3.7 but causes a significant error in r ej 2 .i retrievals. This difference in error 
between the two bands leads to Z\r e , 3 . 7 - 2.1 varying between —5 /iin and 5 /im. It 
is worth mentioning here that the orthogonality of LUT also plays an important 
role in the so-called plane-parallel bias that will be discussed in section 4.4.4. As 
a result, the widths of the r e , 2 .i, r ei 3 . 7 , and the Ar C} 3 . 7 - 2.1 PDFs in Fig. 4.6b are 
quite wide, which suggests that retrieval uncertainties caused by instrument un- 
certainties, ancillary data uncertainties, etc. probably play an important role in 
causing the large variation of Ar e< 3 . 7 - 2.1 for thin clouds in Fig. 4.4. 

As shown in Fig. 4.6, random error in reflectance tends to result in random 
error in r ei 2 .i and r e , 3 . 7 , which in turn leads to a Ar et 3 . 7 - 2.1 with a mean value 
close to zero. Therefore, the random retrieval uncertainties cannot explain why 
r e> 3.7 is systematically smaller than r ei 2 .i, especially when cloud is thick (i.e. high 
signal-to- noise ratio) . 

4.4.2 Vertical cloud structure 

Several studies have shown that ry^.i and r e> 3.7 retrieval have different sensitivities 
to the in-cloud vertical structure of cloud droplet microphysics, leading to spec- 
tral difference in r e retrieval (Platnick, 2000; Kokhanovsky, 2004; Nakajima et ah, 
2010a; Zhang et ah, 2010). The concept of ‘weighting function’ developed by Plat- 
nick ( 2000 ) provides a convenient framework to assess the sensitivity of r e retrieval 
based on the bi-spectral method to the vertical structure of cloud. One of the two 
weighting functions introduced in Platnick (2000) is based on the maximum pene- 
tration depth of photon into cloud. This weighting function, w m (r, r c ), is defined 

1 dR(r) 

R(t c ) dr 


as 


Wm(r,T c ) 


(4.4) 
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Fig. 4.6. An example to illustrate the impact of random error in cloud reflectance on r 
and r e retrievals. In (a) and (b), random error is assigned to both the 0.86-pm band and 
the SWIR band (2.1-pm or 3.7-pm band). In (c) and (d), random error is only assigned to 
the SWIR band. Impacts on r retrievals are shown in (a) and (c). Impacts on r e retrievals 
and Ar e are shown in (b) and (d). See text for details. 


where r c is the total cloud optical thickness, r is the optical depth of a level in 
cloud from cloud top, R{r c ) is the cloud reflectance, and dR(r ) = R(r + dr ) — R(t) 
represents the fraction of all reflected photons that penetrate to a maximum optical 
depth between r and r + dr. As shown in Platnick (2000), the shape of w m (r, r c ) 
can be used to quantitatively interpret the impact of the in-cloud vertical structure 
on MODIS retrievals. Additionally, it is shown that the value predicted from the 
equation 

r* = / r e (r)u; m (r,r c )(ir, (4.5) 

Jo 

where r*, is the expected r e retrieval value based on the vertical weighting, agrees 
to within ±1 pm with the value from full radiative transfer simulation. Note 
that the vertical structure of cloud has little impact on cloud optical thickness 
retrieval. 

Figure 4.7 shows an example of w m (T, r c ) for r e 2 .i (black line) and r e 3.7 (red line) 
for a cloud layer with adiabatic cloud vertical structure, as shown by the dashed 
r e (r) line. The w m (r, r c ) for r e 3.7 peaks at cloud top, while u> m (r, r c ) for r e 2 .i 
peaks lower in the cloud. This difference and Eq. (4.5) suggest that r e 3.7 retrieval 
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Fig. 4.7. Normalized weighting function re m (r, r c ) for r e 2.1 (blue line) and r e 3.7 (red line) 
as a function of cloud optical depth for a cloud with adiabatic vertical structure. Dashed 
line shows the vertical variation of r e as a function of cloud optical depth. 


has more sensitivity to the cloud microphysics at the uppermost part of cloud while 
r e 2.1 retrieval penetrates deeper into the cloud with maximum sensitivity at about 
r ~ 2 or 3. For the case in Fig. 4.7, r e 2.1 and r e 3.7 predicted from w m (r, r c ) based 
on Eq. (4.5) are 15.2 /n n and 16.1 /nn. Both agree well with full retrievals. 

It is important to note that the general shape of the weighting function is not 
sensitive to the detailed cloud vertical structure. Zhang et al. (2010) demonstrated 
that the weighting functions for two cases, one with r e (r ) increasing and the other 
decreasing from cloud top toward cloud base, are quite similar to each other. The 
invariance of the weighting function shape, according to Eq. (4.5), suggests that 
the retrieved r* is largely determined by cloud vertical structure r e (r). And so 
is Ar St 3 . 7 _ 2 .i as far as cloud vertical structure is concerned. As a consequence, 
r e , 3.7 > r e , 2.1 (^•e, 3 . 7 - 2 . 1 > 0) according to Eq. (4.5) when the cloud droplet size 
increases from cloud base toward cloud top, and vice versa when cloud droplet 
size decreases from cloud base toward cloud top. Various microphysical processes 
can affect cloud vertical structure, such as adiabatic growth, entrainment mixing, 
collision-coalescence, and sedimentation. It is therefore important to understand 
the influence of these processes on the vertical structure of cloud microphysics and 
their consequent impact on Ar ej 3 . 7 - 2 . 1 - 

In a classic adiabatic growth model (Brenguier et ah, 2000), cloud r e increases 
monotonically from cloud base to cloud top. In such a case, r ej 3.7 should be larger 
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than r e 2 . i , as shown in the example in Fig. 4.7. Indeed, some recent studies, such 
as Painemal and Zuidema (2011) and Zhang and Platnick (2011), have found 
some positive Ar e , 3 . 7 - 2 . 1 , which is qualitatively consistent with the classic adia- 
batic growth model as shown in Fig. 4.7. However, as mentioned in section 4.3, 
only a small fraction of MODIS MBL cloud observations have positive Ar e , 3 . 7 - 2.1 • 
Moreover, because the adiabatic condensation growth becomes less efficient as the 
droplet grows larger, the droplet growth rate dr e /dz remains relative small over 
the upper portion of the water cloud (Brenguier et ah, 2000). As a result, the mag- 
nitude of Z\r e , 3 . 7 - 2.1 induced by adiabatic cloud vertical structure is expected to 
be small, generally smaller than 2 /mi (Platnick, 2000). 

Cloud top entrainment could lead to a decreasing-with-height r e structure 
at cloud top. Several recent studies speculated that this cloud top entrainment- 
induced structure plays an important role in causing smaller r e , 3.7 retrieval than 
r e, 2 .i (Breon and Doutriaux-Boucher, 2005; Seethala and Horvath, 2010). It should 
be pointed out that the decreasing-with-height r e structure at cloud top assumed 
in these studies seems to indicate a homogenous mixing process (Baker et ah, 
1980). On the contrary, observational studies have actually found more inhomoge- 
neous mixing cases than homogenous mixing cases (Gerber et ah, 2005). A couple 
of recent studies analyzed in situ measurements of MBL clouds off the coast of 
Peru made during the VOCALS-REx (or American Monsoon Systems (VAMOS) 
Ocean-Cloud- Atmosphere-Land Study Regional Experiment) campaign (Painemal 
and Zuidema, 2011; King et ah, 2012). They did not observe the decreasing-with- 
height r e structure at cloud top. In addition, the strong correlation of Ar ej 3 . 7 - 2.1 
with r e> 2.1 and the threshold behavior of Ar e ^. 7 - 2.1 at r ej 2.1 ~ 15 /im (Fig. 4.3) do 
not easily fit into the cloud top entrainment structure argument. 

Warm rain processes (e.g. collision-coalesce) give rise to embryo drizzle drops 
at cloud top, and also make them grow bigger as they fall from the cloud top, 
potentially leading to a decreasing-with-height r e structure from cloud top toward 
cloud base (Berry, 1967; Berry and Reinhardt, 1974; Pruppacher and Klett, 1997). 
Therefore, cloud vertical structure induced by warm rain processes has been argued 
in many studies to be the primary reason for the large negative Z\r e> 3 . 7 - 2.1 (Chang 
and Li, 2002, 2003; Chen et ah, 2007; Nakajima et ah, 2010a, 2010b; Kokhanovsky 
and Rozanov, 2011). Some behaviors of Ar e ^ 3 . 7 - 2.1 shown in section 4.3 seem in 
favor of this argument. For example, it is seen in Fig. 4.3b that Z\r e> 3 . 7 - 2.1 is 
close to zero when r e , 2 .i < 15 /mi and decreases quickly with r e , 2 .i when r e , 2 .i > 
15 /im. This threshold behavior could be explained as a result of increasing drizzle 
probability when cloud r e becomes larger than 15 /mi (Gerber, 1996; Nakajima et 
ah, 2010a, 2010b). Although there is increasing evidence indicating that warm rain 
processes play an important role in causing the large negative Z\r ej 3 . 7 - 2 . 1 , there are 
also studies suggesting otherwise. For example, a couple of case studies based on 
large-eddy simulation of MBL clouds and radiative transfer simulations found only 
a negligible impact of drizzle on both r e retrievals and thereby Ar et 3 . 7 - 2.1 (Zinner 
et ah, 2010; Zhang et ah, 2012). Although it is difficult to draw any statistical 
conclusions from these case studies, they demonstrated a new path toward better 
understanding the impact of drizzle on MODIS r e retrieval that is worthy of further 
exploration. 
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Finally, it should be noted that some /lr e , 3 .r- 2 .i behaviors, the correlation of 
Z\r ej 3 . 7 _ 2.1 with cloud horizontal inhomogeneity index H a shown in section 4.3 in 
particular, cannot be easily explained by cloud vertical structure. This indicates 
that other factors also have a role. 

4.4.3 Cloud droplet size distribution 

As mentioned in section 4.2, in the operational MODIS retrieval algorithm, cloud 
DSD follows the monomodal Gamma distribution with fixed effective variance 
v e = 0.1 (i.e. Eq. (4.2)). Several studies have shown that, when cloud microphysics 
deviates from what is assumed in the MODIS retrieval algorithm, the retrieved r e 
also deviates from the true r e . Invalid microphysics assumptions affect r e , 2 .i and 
7 "e, 3.7 retrievals to a different extent, leading to spectral differences (Chang and Li, 
2001; Minnis et al., 2004; Zhang, 2013). 

The monomodal Gamma distribution assumption for cloud DSD is especially 
problematic when MBL clouds are precipitating, because warm rain processes (e.g. 
collision- coalescence) could broaden cloud DSD and even give rise to a second mode, 
the so-called drizzle or rain mode, creating bi-modal DSD (Berry, 1967; Berry and 
Reinhardt, 1974; Pruppacher and Klett, 1997). In such a case, the MODIS retrieval 
process is to find an r* (the superscript * is to indicate that it is a retrieved value 
under monomodal DSD assumption) that satisfies the following equation: 


where R\ is the observed SWIR band cloud reflectance that is a function of the 
true bi-modal cloud DSD n{r) and Rx^lut^s the SWIR band cloud reflectance 
in the LUT which is pre-computed based on the monomodal DSD assumption. 
Additionally, A denotes the wavelength of the SWIR band used for r e retrieval 
(e.g. 2.1 prn or 3.7 pm). 

Recently, Zhang (2013) developed a semi-analytical model to explain and pre- 
dict MODIS r e retrieval results from Eq. (4.6). In a numerical test, the r e retrievals 
predicted by this model agree well with numerical solutions based on radiative 
transfer simulations (Zhang, 2013). The model is based on the fact that depen- 
dence of R\ on cloud microphysics arises mainly from the dependence of cloud 
single-scattering albedo uj\ on cloud microphysics (Zhang et al., 2009). Figure 4.8a 
shows the 0:2.1 (solid blue line) and 0 ) 3,7 (solid red line) of individual cloud 
droplets computed using a MIE code (Wiscombe, 1979) as a function of droplet 
radius over the range r G [ 1 pm, 1000 pm]. Figure 4.8b shows the so-called bulk 
scattering albedo 0:2.1 ,lut and 0:3.7 Mm ,L{7T hi the MODIS LUT, which are 
averaged over monomodal Gamma DSD, for r e G [5 pm, 30 pm], Zhang (2013) 
expanded 0:2.1 Mm and 0:3.7 Mm into polynomials of r\ 


where ki )X are fitting constants given in Table 4.1, and 0 : 2.1 fim,LUT and 0 : 3.7 ,lut 
into a linear function of r e : 


R\,LUT{r* e ) \ Vb = R\[n(r )] , 


(4.6) 



(4.7) 


(r e ) fc 0 ,A + k l)X r e , 


(4.8) 
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Table 4.1. Fitting coefficients in Eq. (4.7) and the coefficients (values in parentheses) for 
each nonlinear term in Eq. (4.9). The unit of droplet radius in Eq. (4.7) is micron. 


k 0 

ki 

&2 

kj, 

ki 

k 5 

M 

(. ki/ki ) 

( kilkx ) 

{ks/ki) 

(ki/ki) 

(k 5 /ki) 

2.1 pm 9.952E-01 

— 1.736E-03 

2.442E-06 

— 1.187E-09 



(1.17 pm) 

(9.266E-01) 

(— 1.303E-03) 

(6.335E-07) 



3.7 pm 9.623E-01 

— 6.443E-03 

3.348E-05 

— 7.606E-08 

7.757E-11 

— 2.907E-14 

(1.74 pm) 

(9.314E-01) 

(— 4.840E-03) 

(1.099E-05) 

(— 1.121E-08) 

(4.201E-12) 


where fitting constants A; 0 ,a and fci,A are given in Table 4.2. Without going into de- 
tail, Zhang (2013) arrived at the following equation to explain and predict MODIS 
r e retrieval results from Eq. (4.6): 


r 


* 

e,A 


ro, a 


Ah , a 

= T e 

Ah, A 


^ fci.A (^ +2 ) 


(4.9) 


where r* A is the predicted MODIS r e retrieval, (r l ) = J 0 °° r l n{r ) dr is the ith mo- 
ment of the true DSD, and ro,A = {ko,\—ko,\) /k\ is a constant. With the coefficients 
known, r* can be easily solved from the above equation once the ( r 1+2 ) / <V 2 ) terms 
are derived from any given DSD. 


Table 4.2. Fitting coefficients in Eq. (4.8). The unit of droplet radius in Eq. (4.8) is 
micron. 








k 0 

ki 

2.1 

/illl 

9.974E-01 

— 1.874E-03 

3.7 

\i m 

9.744E-01 

— 6.918E-03 


It can be seen from Table 4.1 that ki t \/k\ is close to unity, which suggests 
that r* should be close to the true value of r e , if the higher-order terms, namely 
(Ah, a (r i+2 )) / (k x (r 2 )), in Eq. (4.9) are small. The magnitude of higher-order 
terms is determined by two competing factors. On one hand, the (r l+2 ) / (r 2 ) term, 
which is in the same order of magnitude as r* , increases algorithmically with i. On 
the other hand, as seen in Table 4.1, the Ah,,\/Ah,A term decreases with polynomial 
order i. Provided that r e ranges generally from a few to a few tens of microns, the 
decrease of ki^\/k\^\ terms can be expected to be dominant, leading to a decreasing 
impact of nonlinear terms on r* with increasing order i. Therefore, the most im- 
portant nonlinear term is the second-order term (i — 2). Because k 2 ,\/ki t \ < 0 (see 
Table 4.1) and (r 4 ) / (r 2 ) > 0, the sign of this term is negative, implying that r* 
tends to underestimate the true value when this term is significant. Physically, this 
term is a result of the nonlinear relationship between uJ\ and r. As seen in Fig. 4.8a, 
u\ decreases monotonically with r, and the rate becomes slower as oj\ approaches 
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a) Single-scattering properties vs. r b) Bulk scattering properties vs. r e 



Fig. 4.8. (a) The single-scattering albedo of cloud droplet in the 2.1 - fim (blue) and 
3.7-/im (red) band, and the asymmetry factor in the 0.86-/im band (black) plotted as 
a function of cloud droplet radius. The dashed lines indicate polynomial fitting to the 
single-scattering albedos based on Eq. (4.7). (b) The bulk single-scattering albedo in the 
2.1-/im (blue) and 3.7-/im (red) band, and the asymmetry factor in the 0.86-/rm band 
(black) plotted as a function of cloud droplet effective radius. The dashed lines indicate 
polynomial fitting to the single-scattering albedos based on Eq. (4.8). In the computation, 
monomodal Gamma size distribution with v e =0.1 is assumed. 


the asymptotic value of 0.5. Thus, the curvature of the relationship between oj\ 
and r results in the negative value of the k 2 ,\/ki^\ term. The magnitude of this 
term is determined by the ratio (V 4 ) / (r 2 ) which can be shown to be equivalent 
to r 2 (l + v e ). This relationship reveals that the impact of the second-order term 
depends on r e and to a much lesser extent on v e , namely the width of the cloud 
DSD. This relationship also implies that, in order for the nonlinear term to signif- 
icantly affect UJ\, both r e and v e have to be sufficiently large. Interestingly, such 
conditions happen to be favored in precipitating clouds, as they tend to have larger 
r e and precipitating processes (e.g. collision-coalescence) tend to broaden the DSD 
(i.e. increase v e ) (Berry, 1967; Berry and Reinhardt, 1974; Pruppacher and Klett, 
1997). 

Focusing now on the spectral dependence of Eq. (4.9), it is seen from Ta- 
ble 4.1 that the magnitude of ki^.i^m/ki^.i^m is substantially smaller than 
ki, 3 . 7 fj,m/ki^. 7 fj,m- It indicates a stronger impact of the nonlinear terms on r e re- 
trieval in the 3.7-/im band than the 2.1-/im band. For the retrieval, only the 

second-order nonlinear term is significant. But, for the r* 3 7firn retrieval, higher- 
order terms need also to be included, although the second-order term is still the 
dominant nonlinear term. As discussed above, the second-order nonlinear term 
tends to result in underestimated r* retrieval. This underestimation is expected 
to be stronger in the 3.7-//m band than in the 2.1-/im band, resulting in negative 
r e,3.7 fim ~ r *,2.i jim' This theoretical prediction is consistent with the numerical re- 
sults in Minnis et al. (2004) and Nakajima et al. (2010a). It was found in these stud- 
ies that increasing rain mode in a bi-modal DSD leads to underestimated r e retrieval 
and the underestimation is more severe in the 3.7-//m band than the 2. l-//m band. 

In summary, Zhang (2013) illustrated that r e retrieval-based monomodal DSD 
assumption tends to underestimate the r e of bi-modal DSD, due to the nonlinear 
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relationship between cloud droplet single-scattering albedo and cloud droplet size. 
The underestimation is stronger for r e 3 7(jm than for r e 2 _ i Mm , leading to a negative 

A r e, 3. 7—2.1- 

4.4.4 Plane- parallel r e bias 

As mentioned in section 4.2, one of the fundamental assumptions made in the 
MODIS retrieval algorithm is that a cloudy 1 km x 1 km pixel is horizontally 
homogenous. However, it is known that clouds can have significant horizontal 
variability at smaller scales, leading to sub-pixel inhomogeneity (e.g. Davis et ah, 
1994; Marshak et ah, 1995). When the (horizontally) homogeneous pixel assump- 
tion breaks down, both MODIS r and r e retrievals face challenging issues. For r 
retrieval, for example, sub-pixel inhomogeneity can cause the well-known plane- 
parallel albedo bias. That is, the r retrieval based on the average reflectance of a 
heterogeneous cloudy pixel tends to be smaller than the average of the sub-pixel 
scale r (e.g. Cahalan et ah, 1994; Barker, 1996; Oreopoulos and Davies, 1998). 
Sub-pixel inhomogeneity also affects r e retrieval in various ways, one of which is 
the so-called plane-parallel r e bias (Zhang and Platnick, 2011; Zhang et ah, 2012). 
Plane-parallel r e bias tends to result in overestimation of r e and it affects the 
r e 2 .i retrieval more than r e 3 . 7 , therefore leading to spectral difference (negative 
Zlr e) 3 . 7 _ 2 .i) (Zhang and Platnick, 2011; Zhang et ah, 2012). 

The plane-parallel r e bias is defined as the impact of small-scale variability 
in t on r e retrievals that use area-averaged reflectance (Zhang and Platnick, 2011; 
Zhang et ah, 2012). This bias is illustrated using two idealized examples in Fig. 4.9, 
which shows forward calculations of reflectances in 0. 86 -/ 1111 , 2.1 -/mi and 3.7-/im 
bands. In Fig. 4.9a and 4.9b, we assume that half of a MODIS pixel overlying a 
black surface is covered by a cloud with n = 2.8 and r e — 8 //m and the other 
half is covered by a cloud with T 2 = 30.8 and r e — 8 //in. Panels c and d assume 
the same optical thickness but use r e — 18 /mi. Focusing on the r retrieval, the 
figure illustrates the well-known ‘plane-parallel albedo bias’ (Cahalan et ah, 1994): 
the retrieved r based on the mean reflectance of inhomogeneous pixels tends to be 
smaller than the mean of the sub-pixel r. In this example the value of retrieved 
from the averaged reflectance (r = 10 . 8 ) is substantially smaller than the average 
value (r = 16.8). 

This problem is more acute for retrievals of r e because the reflectance used to 
infer r e also depends on r over much of the range of plausible values, as has been 
shown in section 4.4.1. If the reflectance at non-absorbing and absorbing wave- 
lengths depended only on r and r e , respectively (i.e. if the LUT was orthogonal), 
reflectance at absorbing wavelengths would be uniform in our example and parti- 
cle size could be retrieved perfectly. As Fig. 4.9 demonstrates, however, the LUT 
is not orthogonal. The nonlinearity leads to a simultaneous underestimation of r 
(plane-parallel-albedo bias) and overestimation of r e (plane-parallel r e bias). The 
area over which this is true is larger in the less-absorbing band, which explains why 
r c , 3.7 is smaller than /y. 2 , 1 . The impact becomes more pronounced as r e increases: 
in Fig. 4.9c and 4.9d, the true r e — 18 //m while r e , 2 .i and r e> 3.7 retrieved from 
averaged reflectances are 24 /mi and 18.1 /mi, respectively, resulting in a Ar e ^. 7 - 2.1 
around —6 /mi. 
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Fig. 4 . 9 . Two theoretical cases to illustrate the nonlinearity effect in r e retrievals resulting 
from sub-pixel cloud inhomogeneity. Numbers on top of the Nakajima-King look-up-table 
(LUT) curves correspond to values of r contour lines in the LUT and the numbers on the 
right of the curves correspond to values of r e contour lines in the LUT. 


Recently, Zhang et al. (2012) developed a numerical MODIS retrieval simulator 
based on the combination of an LES model with bin microphysics (Ackerman et ah, 
2004; Fridlind and Ackerman, 2010) and radiative transfer models (both ID and 
3 D). In the case study based on this simulator, they confirmed that the r e plane- 
parallel bias, as conceptually illustrated in Fig. 4.9, indeed play an important role 
in causing the spectral difference between r 6) 2.1 and r e , 3.7, especially for those pixels 
with high sub-pixel inhomogeneity. 

As shown earlier in Fig. 4.5a, the statistics of MODIS r e> 2.1 show a strong 
dependence on the sub-pixel inhomogeneity index, H a . When a cloud is relatively 
homogeneous (i.e. H a < 0.3) the peak of the PDF of MODIS r e , 2 .i is around 
12 ~ 14 pm. When H a becomes larger than about 0.3, the peak increases quickly 
with H a , by up to about 20 pm when H a is about unity. Interestingly, MODIS r e , 3.7 
shows little dependence in Fig. 4.5b, leading to significant negative Ar Cj 3. 7-2.1 when 
H a > 0.3. The behaviors of r e , 2 .i> r e,3.7 and Ar e , 3 . 7 - 2.1 in Fig. 4.5 are consistent 
with the concept of r e plane-parallel bias as shown Fig. 4.9, therefore suggesting 
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r e plane-parallel bias is an important factor causing the Z\r e , 3 . 7 - 2 . 1 - On the other 
hand, Fig. 4.5 also suggests that the impact of r e plane-parallel bias is limited to 
those highly inhomogeneous pixels, with H a larger than about 0.3. 


4.4.5 3D radiative transfer effect 

As shown in a couple of recent studies (Zhang and Platnick, 2011; Zhang et ah, 
2012), 3D radiative transfer effects generally influence r 6 j 2.1 and r 6j3 7 retrievals to 
different extents. This makes 3D radiative transfer a potentially important factor in 
causing the spectral dependence of MODIS r e retrieval as described in section 4.3. 

As mentioned in section 4.2, the MODIS retrieval algorithm assumes the so- 
called plane-parallel cloud model. However, clouds in reality have significant hor- 
izontal structures and variabilities at various scales, from thousands of kilometers 
to a few meters (Cahalan et ah, 1994; Davis et ah, 1994; Wood and Hartmann, 
2006; Wood et al., 2008). Horizontal cloud variability gives rise to horizontal photon 
transport (Davis and Marshak, 2010), which makes the observed reflectance depen- 
dent on not only the property of pixels within the field of view, but also the proper- 
ties of surroundings (Marshak et ah, 1995). In such circumstances, cloud properties 
retrieved under the plane-parallel cloud assumption are subject to significant er- 
rors. This is the so-called 3D radiative transfer effect. Take the so-called ‘shadowing 
effect’, for example (Varnai and Marshak, 2002; Marshak et ah, 2006). Consider, 
for example, a cloudy pixel that is in the shadow of an adjacent pixel. The dark- 
ening of the pixel would be fallaciously interpreted under the plane-parallel cloud 
assumption, leading to underestimated r and overestimated r e retrievals (Varnai 
and Marshak, 2002; Marshak et ah, 2006; Zhang and Platnick, 2011; Zhang et ah, 
2012 ). 

It was shown in a couple of recent studies (Zhang and Platnick, 2011; Zhang et 
ah, 2012) that 3D radiative transfer effects generally have a stronger impact on the 
r c , 2 . 1 retrieval than on the r e> 3.7 retrieval. For example, when the shadowing effect 
leads to an overestimation of r e , the magnitude of overestimation is larger for r e , 2 .i 
than for r e , 3.7. The reason for this is two-fold. First, 3D radiative transfer effect 
results from horizontal photon transport and the strength of droplet absorption in 
the 3.7-/im band exerts a stronger limit on horizontal photon transport than the 
weaker 2 . 1 -/im band absorption. Secondly, as shown in sections 4.4.1 and 4.4.4, 
the LUT for the 0.86-//m and 2. l-//m combination is less orthogonal than that for 
the O. 86-//111 and 3.7-//m combination. As a consequence, the r e , 2 .i retrieval is more 
sensitive to 3D radiative transfer effects on the 0. 86-^1 band reflectance than r e , 3.7 
retrieval. 

Recently, Zhang et al. (2012) studied the impact of 3D radiative effects, as well 
as other factors, on MODIS cloud property retrievals using an LES model with 
a binned microphysical scheme (Ackerman et ah, 2004; Fridlind and Ackerman, 
2010) coupled with a 3D radiative transfer model (Cahalan et ah, 2005; Pincus and 
Evans, 2009). Figure 4.10 shows one of the LES cases in this study, which is based 
on an idealized case study (Stevens et ah, 2001) from the Atlantic Trade Wind 
Experiment (ATEX), with an average cloud droplet concentration (weighted by 
liquid water mixing ratio) of about 30 cnH 3 . For more details about this LES case 
and the configurations of radiative transfer models, readers are referred to Zhang 
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et al. (2012). Figure 4.10a provides a planar view of the cloud r derived from the 
LES for this case. Figure 4.10b shows the cloud extinction coefficient along the 
vertical cross-section indicated by the red dashed line in Fig. 4.10a. In radiative 
transfer simulations, solar zenith angle is assumed to be 20° and viewing zenith 
angle to be 0°. The ID and 3D simulations of cloud reflectances along the cross- 
section for the 0.86-/im, 2.1-zim and 3.7-yum bands are shown in Figs. 4.10c, 4.10d 
and 4.10e, respectively. Although ID and 3D simulations are in general agreement, 
the influence of horizontal photon transport can be clearly seen in a few regions, 
especially in Fig. 4.10c. For example, at x — 7 km the 0.86-/im cloud reflectance 
based on ID simulation is substantially larger than 3D simulation. But, on both 
sides of this reflectance peak (i.e. x — 6.5 km and x — 7.5 km), 3D reflectance is 
larger than ID. These differences between ID and 3D simulations are the result of 
horizontal photon transport due to multiple scattering (Varnai and Davies, 1999). 
Interestingly, this horizontal photon transport effect is much reduced in the 2.1-/im 
and 3.7-//ni simulations. This is because the absorption in the SWIR bands limits 
multiple scattering and therefore reduces the horizontal transport of photons. 

Figure 4.11 shows the impacts of 3D radiative transfer effects on r e , 2.1 and 
r e> 3.7 for the case in Fig. 4.10. In the figure, we use the difference between the 
r retrieval based on 3D radiative transfer simulation (referred to as ‘3D r’) and 
that based on ID radiative transfer simulation (referred to as ‘ID r’) as an index 
of the 3D radiative transfer effects. By definition, 3D r > ID r in the case of 
illumination, and 3D r < ID r in the case of shadowing. As expected, the r e 
retrievals based on 3D radiative transfer simulation (‘3D r e ’) appear smaller than 
those based on ID radiative transfer simulation (‘ID r e ’) in the case of illuminating 
effect (i.e. 3D r e — ID r e < 0), and vice versa in the case of shadowing effect (i.e. 
3D r e — ID r e > 0). What is interesting in Fig. 4.11 is that the impact of 3D radiative 
effects on r e ^.i is generally stronger than on r e> 3.7, for the reasons discussed earlier 
in this section. Another important point to note in Fig. 4.11 is that illuminating 
and shadowing effects naturally come in pairs and tend to cancel each other out 
(Marshak et ah, 2006). Therefore, 3D radiative effects seem unlikely to cause large 
systematic bias in r e retrieval. Moreover, since 3D radiative effects are a result of 
horizontal photon transport, its impact on MODIS retrievals is limited to pixels 
with large horizontal inhomogeneity. In addition, when a cloud pixel has strong 
horizontal inhomogeneity, both 3D radiative effects and the plane-parallel r e bias 
discussed in section 4.4.4 may play a role at the same time, creating either positive 
or negative ^, 37 - 2.1 (Zhang et ah, 2012). 


4.5 Discussion 

In the previous section, we reviewed several mechanisms that could cause significant 
difference between MODIS r e , 2.1 and r e; 3.7 retrievals for MBL clouds. But the most 
fundamental reason is that MBL clouds in reality do not always satisfy the ideal 
cloud model (e.g. plane-parallel and homogenous) on which MODIS is based. When 
the cloud deviates from the ideal model, the deviation affects the r e , 2 .i and r e , 3.7 
retrievals to different extents or in different ways, leading to significant difference 
between the two. From this perspective, spectral dependence is an inherent and 
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Fig. 4.10. (a) Planar view of the cloud r of the ATEX clean case at 6 h of simulation 
time, (b) The cross-section of cloud extinction coefficient (/3) along y — 2 km in (a). Cloud 
bi-directional reflectance along the cross section is shown for the (c) 0.86-/xm, (d) 2.1-fim, 
and (e) 3.7-y.m MODIS bands simulated using ID (blue) and 3D (red) radiative transfer 
models. 
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inevitable feature of the r e retrieval based on the bi-spectral method. On one 
hand, it reflects the limitations of the method but, on the other, it indicates that the 
operational MODIS retrieval results are generally consistent with our expectations. 

4.5.1 Which one is better? 

If spectral dependence is inevitable in MODIS r e retrieval, then a natural question 
arises: which one, r e 2.1 or r ej 3.7, is better or more correct? Unfortunately, there is 
no simple answer to this question. As shown in section 4.4.1, r ej 3.7 is less affected 
by random retrieval uncertainties, if uncertainties associated with 2.1-/im and 3.7- 
pm band observations are of similar magnitude. Also, in comparison with r e 2.1 , 
r e,3.7 is generally less affected by the plane-parallel r e and 3D radiative effects as 
discussed in sections 4.4.4 and 4.4.5. In these cases, r e< 3.7 is expected to perform 
better than r 6 ) 2 . 1 - However, as discussed in section 4.4.2, from the perspective 
of vertical weighting, 2.1 and r 6) 3.7 simply reflect different parts of the cloud. 
Simply put, there is no right or wrong. Moreover, as discussed in section 4.4.3, 
when the true cloud DSD is bi-modal, both r G) 2.1 and r ej 3.7 tend to underestimate 
the true r e , but the underestimation is more severe in the 3.7 than in the r e> 2 . 1 - 
Therefore, whether to use r e> 2.1 or r e> 3.7 depends on the intended application and 
the geographical location of interest. For example, several studies have attempted 
to retrieve the vertical profile of MBL clouds utilizing the combination of r (: 2.i 
and r e , 3.7 retrievals from MODIS (Chang and Li, 2002, 2003; Kokhanovsky and 
Rozanov, 2011). For these applications, both r e , 2.1 and r ei 3.7 are needed and their 
difference serves as a useful signal rather than a data issue. As shown in Fig. 4.2, 
r 6 j 2 .i and r 6j 3.7 agree quite well over the coastal stratocumulus regions. So, if those 
regions are of interest, it perhaps makes little difference which retrieval is used. 
In contrast, over broken cloud regions where clouds have significant horizontal 
heterogeneity, r e , 3.7 would be less affected by plane-parallel r e bias and 3D radiative 
effect and therefore might be a better choice than r ej 2 . 1 - 
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It should be noted here that r e> 3.7 is reported as having a respective difference 
to successful r e , 2.1 retrieval in the Collection 5 of MODIS cloud product. Thus, in 
Collection 5, the sampling of r e , 3.7 is biased by the success of r e , 2 . 1 - In addition, 
although retrieved in Level 2 products, r ej 3.7 is not aggregated to Level 3 products. 
In the current Collection 6 of MODIS cloud product, these sampling biases will be 
removed and r e , 3.7 will be sampled independently and aggregated to Level 3 data. 
These changes will facilitate the use of r e , 3 . 7 . 

4.5.2 Cloud regime classification 

Summarizing the discussions in the previous section, one can see that none of 
the hypotheses can uniquely explain all of the Ar S} 3 . 7 - 2.1 features shown in sec- 
tion 4.3. This is surprising because all these mechanisms are more or less entangled 
together in the retrieval. For example, when a cloud pixel has strong horizontal 
inhomogeneity, both the plane-parallel r e bias (section 4.4.4) and 3D radiative ef- 
fects (section 4.4.5) may play a role at the same time. When warm rain processes 
begin to develop in MBL clouds, it can lead to both cloud vertical structure and 
bi-model DSD, in which case both the vertical weighting effect (section 4.4.2) and 
the DSD sensitivity effect (section 4.4.3) will influence MODIS r e retrievals. Thus, 
for a better understanding, we need to untangle these mechanisms and sort out 
their relative importance in different cloud regimes. By ‘cloud regimes’, we mean 
groups of MBL clouds with different characteristics in terms of optical (e.g. r), mi- 
crophysical (e.g. r e ), and macrophysical (e.g. horizontal properties). An attempt of 
such cloud regime classification and analysis has been made in Fig. 4.12. It shows a 
color contour of the mean value of Ar e ^. 7 - 2.1 projected on sub-pixel heterogeneity 
index H a and r 6 j 2.1 for MBL clouds with r > 5. The black contour lines indicate 
the joint frequency histogram of H a and r ei 2 .i based on one month of MODIS ob- 
servation. Based on the combination of r e 2 .i and H a , MBL cloud pixels observed 
by MODIS can be classified into the following three regimes: 

Regime 1 : MBL cloud pixels in this regime have low horizontal heterogeneity 
(H a < 0.3) and relatively small effective radius (r 6 i 2.1 < 20 pm). Low H a 
index suggests that the plane-parallel r e bias and 3D radiative effects should 
be small in this regime. Small r e> 2.1 suggests that the probability of warm rain 
precipitation is relatively small because the critical value of r e for collision- 
coalescence to occur is around 15 pm (Gerber, 1996; Rosenfeld et ah, 2012). 
As a result, the vertical weighting effect and the DSD sensitivity effect can be 
expected to be small. It now becomes clear that MBL cloud pixels in this regime, 
with low horizontal heterogeneity and low possibility of precipitation, satisfy 
those assumptions made in MODIS cloud retrieval mentioned in section 4.2. 
For this reason, it is not surprising to see good agreement between r ej 2 .i and 
r e , 3.7 in this regime (Z\r e; 3 . 7 - 2.1 is within ±2 pm), although a small difference, 
probably due to random retrieval uncertainties, may still exist. 

Regime 2: MBL cloud pixels in this regime are characteristic with high hori- 
zontal heterogeneity ( H a > 0.3). Retrieval uncertainties, plane-parallel r e bias, 
and 3D radiative effects may all play a role in this regime. Note that these ef- 
fects have different impacts on r e retrievals and therefore Ar ej 3 . 7 - 2 . 1 - Retrieval 
uncertainties, as discussed in section 4.4.1, are most likely to cause random 
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Fig. 4.12. The color contour of monthly mean Z\r e , 3. 7-2.1 for clouds with r > 5 on the 
space specified by r e , 2.1 and H a . The black lines indicate the relative frequency of each 
grid box, specified by certain combinations of r e , 2.1 and H a (unity corresponds to the 
most frequently observed combination of r e , 2.1 and Ha-). See text for details on cloud 
regime classification and implications for MODIS effective radius retrievals. 


errors, rather than bias, in Ar e) 3 . 7 - 2 . 1 - Similarly, as discussed in section 4.4.5, 
3D radiative effects can cause either positive or negative Ar e , 3 . 7 - 2 . 1 , depend- 
ing on the nature (e.g. illuminating or shadow) of the 3D effect. In contrast, 
the plane-parallel r e bias, as discussed in section 4.4.4, tends to cause negative 

Ar e , 3 . 7 - 2 . 1 - 

Regime 3: MBL cloud pixels in this regime have low horizontal heterogeneity 
{H a < 0.3) and relatively large effective radius (r ej2 . 1 > 20 ji m). Large effective 
radius suggests that MBL clouds in this regime are possibly precipitating. The 
warm rain process can give rise to the vertical weighting effect and the DSD 
sensitivity effect, as discussed in section 4.4.2 and 4.4.3 respectively. Because 
both effects tend to make r e! 2 .i larger than r 6 j 3 . 7 , thus the bias of this regime 
is large and negative, as seen in Fig. 4.12. 


4.6 Outlook of future work 

Recently, the spectral dependence of MODIS r e retrievals for MBL clouds has re- 
ceived significant and increasing attention. A number of recent studies reviewed 
in this chapter have significantly improved our understanding of the nature and 
potential causes of this issue. Nevertheless, many outstanding questions remain 
and future work is needed particularly in the following areas: (i) More quantitative 
understanding : Together , those hypotheses reviewed in section 4.4 provide a quali- 
tative explanation for the r e , 2 .i and r ej 3.7 differences. However, a more quantitative 
understanding is still lacking. It remains unclear whether the threshold behavior 
of Ar e , 3 . 7 _ 2 .i at r e , 2 .i ~ 15 /im in Fig. 4.3 is a coincidence or a result of warm rain 
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process. In order to answer questions like this, hypotheses must be evaluated in a 
more quantitative way in future work, (ii) Independent measurements : For more 
objective and quantitative evaluation of MODIS r e , 2.1 and f'eM .7 retrievals, inde- 
pendent measurements are needed. For example, several studies have investigated 
the impact of 3D radiative effects on cloud masking (Zhao and Di Girolamo, 2006), 
cloud optical thickness, and cloud effective radius retrievals (Marshak et ah, 2006) 
using the high spatial resolution ASTER (Advanced Spaceborne Thermal Emis- 
sion and Reflection Radiometer) observations. The ASTER observation may prove 
useful in future analyses to understand the relative importance of the impacts of 
3D radiative effects and plane-parallel r e retrieval bias on cloud property retrievals. 
Several studies have demonstrated that high temporal-spatial resolution cloud DSD 
measurements from air-borne in situ instruments, collocated with MODIS obser- 
vation, are very helpful for understanding the errors and uncertainties in MODIS 
retrievals. Other remote sensing measurements, in particular those from A- Train 
satellite sensors (e.g. CloudSat, POLDER, and CALIPSO), would add unique and 
valuable perspectives on MBL cloud microphysics that are worthy of exploration 
in future work, (iii) Combining observations with numerical models : Every re- 
mote sensing method has its limitations. Even though a combination of various 
methods, such as those available from A-Train sensors, provides a complementary 
perspective, there will always be gaps that cannot be covered by observations. As 
shown in a number of recent studies, a combination of LES models with satellite 
retrieval simulators provides a flexible and powerful tool for understanding the 
impact of various factors on MODIS cloud property retrievals (e.g. Zinner et ah, 
2010; Zhang et ah, 2012). This new avenue of research should be further explored 
in future work, (iv) Implication studies : MODIS cloud products are popularly used 
in climate change studies, aerosol indirect effect studies, and climate model valida- 
tions. The implications of the spectral dependence of MODIS r e retrievals for these 
‘down-stream’ applications should be investigated in future work, (v) Novel use of 
r e, 2 .i and r e , 3.7: As discussed in sections 4.4.2 and 4.4.3, in certain circumstances, 
the difference between r e , 2.1 and r ej 3.7 actually contains useful information about 
MBL clouds. Actually, a number of studies have attempted to retrieve the vertical 
profile of precipitating MBL clouds utilizing the combination of effective radius 
retrievals from different spectral bands (Chang and Li, 2002, 2003; Kokhanovsky 
and Rozanov, 2011). Further work is needed along these lines to fully explore the 
potential of MODIS observations. 
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